April 12, 2010 0:34 WSPC/Trim Size: 9in x 6in for Proceedings gerrard3v2 



o 

(N 

Oh 

< 
0\ 



$-1 



X 



EARTH FLYBY AND PIONEER ANOMALIES 



M.B. GERRARD* AND T.J. SUMNER 

Department of Physics, Imperial College London, Prince Consort Road, London. 
SW7 2BW, UK. E-mail: t.sumner@imperial.ac.uk 



o 

n^ Applying Newtonian dynamics in five dimensions rather than four, to a universe 

that is closed, isotropic and expanding, suggests that under certain circumstances 

an additional and previously unidentified acceleration, a a, can arise affecting the 

four dimensional motion of spacecraft. The two cases of this acceleration being 

either real or virtual are considered. In the real case, simple estimates of aj\ 

^f^ are shown to be in partial agreement with reported acceleration anomalies from 

^ several Earth flybys and from the Pioneer spacecraft. However, these estimates do 

QQ not fully reconcile with radio Doppler tracking data. The virtual case, by contrast, 

iy~\ appears to overcome these and other difficulties with the real case, and is discussed 

y—{ in an addendum. Furthermore, the virtual case has an altitude dependence which 

rf^ makes detection of any anomaly unlikely above '~2000km. Equations governing 

• this additional acceleration have been derived from first principles, without the 

|- _ introduction of free parameters or new constants and without amendment to the 

^-^ law of gravity. 

00 

o 
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J, Measurements of the velocity profiles of several spacecraft during Earth fly- 

C^ bys, and of Pioneers 10 and 11 following their flybys of, respectively, Jupiter 

and Saturn have identified small discrepancies between expected and ob- 
served velocities [1,2]. In particular for Earth flybys, relative discrepancies 
in velocity of the order 10^® (which can be either positive or negative) have 
been reported and, for the two Pioneer spacecraft, a near constant accel- 
eration towards the Sun of the order 8.74 x 10""'^° ms~^ has been detected. 
Currently, these anomalies are without widely accepted explanations. 

The purpose of this paper is to show that the introduction of a fifth 
large-scale dimension to the Newtonian equations of motion may simulta- 
neously account for both these anomalies. 



1. Introduction 
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2. A Five Dimensional Framework 

An additional large-scale dimension, r, that is orthogonal to the three space 
dimensions s{x, y, z) and to the time dimension, t, is introduced and iden- 
tified as the radius of curvature of four dimensional space-time in a closed, 
isotropic and expanding universe. The dimension r is defined to have units 
of metres and a background value remote from gravitating bodies of r„. 
The numerical value of r„ is not required for the current analysis, although 
one has been proposed by the authors [3]. Radiation freely propagating 
at the speed of light, c, in this closed isotropic universe will experience an 
acceleration, a^, in the direction of r and varying as a^ — <? jr, for a five 
dimensional universe that is locally fiat (sec Fig. 1)^. To assist visualisa- 
tion of the relationship between space, time and the additional dimension, 
r, only a single space dimension, s, is shown in Fig. 1. Using spherical 
polar co-ordinates, the relationships between s, t and r are given by: 

Ts = rsina (1) 

s = 0r, (2) 

cdt — rda (3) 

From these three equations it is possible to derive an expression for the 
background expansion of an arc of space, s, as follows: 

ds c ,,. 

^s-cota (4) 

at r 

For each point in space, an angular velocity vector Via can be defined 
which subtends an angle a relative to the dimension r and whose value is 
given by equation (5) in which c is defined to be positive: 

l^s = - (5) 

The presence of a gravitating body of mass M causes the radius of cur- 
vature of space-time, r, to reduce relative to its limiting value r„, according 
to the following approximation [3] : 



GM \ 



J (6) 



^Similarities have been identified by the authors between ar and the MOND parameter 
ao [3,4] 
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Figure 1. An arc of space, s, in a closed, isotropic and expanding universe 



The {s,r) contour defined by this equation is the boundary between four 
dimensional space-time and the fifth dimension r. It is the locus of points 
along which the resolved accelerations of gravity gg (towards M) and of a^ 
(in the direction of r increasing) are equal and opposite. Hence for r c^i r^'- 

dr GM 



9s 



ds 



(7) 




Figure 2. The influence of a gravitating body, M, on the radius of curvature of space- 
time, r, and motion of a small body, m, along the associated contour r (M, s) 



For a small body (shown as m in Fig. 2) travelling at a relative radial 
velocity, Vradj with respect to M, the derivative of the angular velocity 
vector fig is given by equation (8). At this point in our analysis only 
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variations in r^ due to the motion of the body m through the gravitational 
field Qg are considered (i.e. a is assumed constant, so ignoring the expansion 
of the universe): 



1 V. 



Vis = —gs (8) 



c 



Clg gives rise to an additional acceleration a a in the space dimension, s, 
as defined by the vector cross product in equation (9): 

aA^Ts A Vis (9) 

Substituting from equation (8) provides the following expression for the 
acceleration a^: 

aA^gs (10) 

c 

a A acts along the line between M and m, i.e. is co-aligned with gg. It 
may alternatively be expressed, by substituting from equation (7), solely in 
terms of the fifth dimensional parameter, r, as shown in equation (11) from 
which it can be seen that qa is a five dimensional analogue to the Coriolis 
acceleration of four dimensional dynamics. 

aA = --Ar (11) 

r 

The additional acceleration, a a, is not a modification to the gravitational 
vector 5s, so any change in the kinetic energy of the mass m under the 
action of qa is not offset by a reduction or gain in gravitational potential 
energy of m with respect to M. Accordingly, the conservation of energy 
requires that the action of a^ will only give rise to an observed change in 
the kinetic energy of m where energy may be transferred to/from m by 
virtue of the general dynamics of the system within which it is moving, for 
example in the circumstances of fiybys. 



3. Earth Fiybys 

It is straightforward to show from equation (10) that where energy transfer 
to/from m can take place, the additional change in velocity due to ga 
as m moves radially with respect to M between si and S2, is given by 
equation (12), where 9 is the mean value (between si and S2) of the angle 
subtended by the radial velocity vector, Vrad, to the direction of motion of 
the spacecraft, such that Av = Avrad cos 9 - noting that the acceleration 
a A and the radial vector are always aligned (see Fig. 3): 
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GM ( \ 1 \ - 

Aw= COS0 (12) 

C Vsi 52/ 

However, in the case of planetary flybys, it is the vector cross-product 
between the angular momentum of the planet, in this case the Earth, around 
the Sun and the motion of the spacecraft which defines the energy transfer. 
Accordingly, the velocity vector of the spacecraft (w) must be resolved into 
the plane of the ecliptic and it follows that the cosine of the angle /?' - the 
ecliptic angle of incidence of the velocity vector (v) - determines the extent 
to which energy transfer takes place as a consequence of a a acting on the 
spacecraft. Similarly, the angle of divergence d' between the velocity {v) 
and radial velocity (v^ad) vectors projected into the ecliptic (rather than 
simply ff) determines the relationship between a a and changes in velocity 
of the spacecraft, such that Aw oc a^ cos /?' cos 0' . The change in velocity 
(Aw) is then given by equation (13) where Mg is the mass of the Earth: 

Aw = ^^ (- - 1^ cos/3'cosF (13) 

C V*l *2/ 

For an inbound trajectory towards Earth, Aw^ will be positive and, 
conversely, for an outbound trajectory Awq will be negative. Hence, the net 
gain/loss in kinetic energy of the spacecraft during a flyby due to a a will 
be given by the aggregate of inbound and outbound differences in velocity 
(i.e. Awi + Awo). For a flyby measured over a distance from/to s via the 
distance of closest approach (given the symbol sy), the additional velocity 
gain/loss on flyby is given by: 

Aw = Aw, + Awo = \ I (cos/3'cos0'j - cos /3' cos F^) (14) 

c \si s) 

There is an important short-coming to this expression in that during 
the flyby the angle of incidence /?' of the velocity vector to the ecliptic is 
not constant, particularly during the transition from inbound to outbound 
trajectories. This will tend to introduce significant inaccuracy for trajec- 
tories that are highly asymmetrical about the plane of the ecliptic. By the 
same token, use of an average value, 6'', rather than the full time-dependent 
profile Q' [€) during the flyby will also be a source of inaccuracy. The char- 
acteristics of equation (14) arc such that for the six reported Earth fiybys 
(Galileo(1990), Galileo(1992), NEAR(1998), Cassini(1999), Rosetta(2005) 
and MESSENGER(2005)), the key portion of the trajectories which gives 
rise to the net additional velocity gain/loss Aw is roughly 2 hours either 
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Figure 3. Divergence between velocity, v, and radial velocity, v,.ad vectors during flybys 



side of the closest approach. This arises for two reasons: first, the mag- 
nitude of aA falls away with distance from Earth (s) and with a declining 
radial velocity {vrad) of the spacecraft, such that for each of the inbound 
and outbound trajectories, the change in velocity (Aw) over the 2 hour 
segment measured from the closest approach, is circa 90% of that arising 
over the respective 4 hour interval (and is circa 80% of that arising over 
the respective 24 hour interval); and second, a a influences the velocity of 
the spacecraft (u) via the assymetric time dependent profile of the vector 
cross-product of the flyby itself. 

To take proper account of the time dependency of /3' (t) and 9' (t) re- 
quires detailed numerical integration and modelling of the trajectories dur- 
ing the sectors -I-/- 2 hours either side of the closest approach, which is 
beyond the scope of this paper. However, a simple approximation can be 
introduced by splitting the inbound and outbound calculation into just two 
sectors each - the first is from the 2hr boundary location, S2Hr, in to 2s f 
(the outer sector), and the second for distances Sf < s < 2s f (the inner 
sector). In the outer sector, 9 and so 9' is assumed zero for all inbound and 
outbound trajectories and in the inner sector during which 30° < 9' < 90°, 
9' is set to an average value of 60° except where a trajectory has a signifi- 
cant value of P' . In these cases only (NEAR outbound and MESSENGER 
inbound) a reduction in 9' needs to be made and a value of 30° has been 
assumed. Also for these two trajectories (NEAR and MESSENGER) j3' 
within the inner sector is given a value half that of the outer sector, as ap- 
propriate, to allow for the highly asymmetric transition of /3' at the closest 
approach in these cases. 
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Table I: Comparison of predicted and observed values for addi- 
tional velocity gained/lost on flyby 



Parameter 


Units 


Galileo 


Galileo 


NEAR 


Cassini 


Rosetta 


M'GER 


Date 




1990 


1992 


1998 


1999 


2005 


2005 


^perigee 


km 


960 


303 


539 


1175 


1956 


2347 


s/ 


km 


7331 


6674 


6910 


7546 


8327 


8718 


S2Hr 


kmxlO-5 


0.73 


0.73 


0.61 


1.19 


0.46 


0.47 


Vf 


kms^^ 


13.740 


14.080 


12.739 


19.026 


10.517 


10.389 


I3i 


Degree 


10.89 


-17.89 


2.44 


-2.12 


2.89 


52.37 


fi'o 


Degree 


-17.66 


-6.71 


-59.86 


-1.62 


-12.43 


-13.73 


Av 


mms^^ 


3 


-5 


13 


-0.03 


2 


[5] 


AVobs 


mms^^ 


3.92 


-4.6 


13.46 


<-0.2 


1.8 


0.02 



As can be seen from Table I, ^ equation (14) applied for simplified two- 
sector trajectories predicts values for Av which are in reasonable agree- 
ment (both in regards to magnitude and sign) with most observations for 
the spacecraft. The degree of agreement can be further illustrated by solv- 
ing for a common distance (s) which, if applied in equation (14) instead 
of the common time-cutoff (±2hrs either side of the perigee) would limit 
the differences between Av and Avobs, particularly for MESSENGER. This 
common distance (s) lies in the range 72 — 76 x 10'^ km and provides agree- 
ment between Aw and Avobs to an accuracy of ±1 mm/s for all spacecrafts^ . 

Alternatively, if during the inner sector for all six trajectories, account 



Values for: (i) the altitude of closest approach (-ffperigee); (ii) the velocity at closest 
approach (vf); and (iii) the measured values for the additional velocity (Avabs) a-fc all 
taken from [1] with the exception of the reported value of Av^^b^ for Cassini which is taken 
from [2], Other values in the table are derived as follows: (a) Sf is the sum of Hpgrigee 
and the Earth's mean radius REarth (6371 km); (b) the distance of the spacecraft from 
Earth approximately 2hr either side of the flyby (s2//r) is taken from [5]; (c) the ecliptic 
angle of incidence (/3') of the velocity vector is approximated by the ecliptic latitude 
(/3) of the spacecraft where the angle of divergence (8) between the velocity (v) and 
radial velocity (v^ad) vectors is close to zero (i.e. other than during the close flyby 
itself), where (/3) is calculated from the equatorial latitude, the right ascension [1] and 
the obliquity of the ecliptic (23.45°); (d) Av is calculated from equation (14) using a two 
sector approximation as explained in the main text. 

'^In fact, the functional form of equation (14) may give rise to better fits of the Doppler 
tracking data to the equations of 4-dimensional dynamics than expected, notwithstand- 
ing the action of a a, for the same reason that equation (22) does in the case of a a being 
interpreted as virtual (see section 6.3). 
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is taken of the change in (3' at the closest approach as the spacecraft tran- 
sition from inbound to outbound trajectories, typically by assuming, where 
appropriate, that the average value for /3' is reduced by 50% then close 
agreement is achieved between observed and predicted values for Aw with 
values of 0' within 60° ± 2° in the inner sector, except in the case of NEAR 
and MESSENGER for which Wo and W^ respectively are reduced to 30° ±5° 
to account for the significant ecliptic angles of incidence of these trajectories 
as noted earlier. 

4. Pioneer Anomaly 

Returning to equation (5) to apply this to the Pioneer Anomaly, we need 
to generate the time derivative of the angular velocity vector, i7s, including 
not only the effect of a body moving through a gravitational field {gs), but 
also the variation of r^ (and so of f2s) with time, due to the background 
expansion of the universe - i.e. a is no longer constant. In the context of 
the Pioneer Anomaly g^ is the gravitational acceleration due to the Sun. 
The time derivative of ils is then given by equation (15): 

• 1 Vrad 1 c2 

lis — — Qs cot q; (15) 

From equation (4), Hubble's constant {K) is identified as the term 
(-cotof). Substituting this and fi,, from equation (15) in equation (9) 
provides a general result for the additional acceleration a^ as shown in 
equation (16). Strictly speaking, the fifth dimensional parameter, r, (i.e. 
the radius of curvature of spacetime) is also itself expanding as a result of 
the conservation of energy at the boundary of four dimensional spacetime'*, 
and in principle this gives rise to a third term in equation (15). How- 
ever, the effect is small compared with the other terms of equation (15). 
In the context of Earth flybys, the second term of equation (15) is itself 
much smaller than the first, which justifies the use of only the first term in 
equation (10). 

aA^gs-^ cB (16) 

c 

The Pioneer Anomaly has been seen as the spacecrafts are travelling 

away from the Sun, so gs is negative with respect to the positive radial 

velocity, v^-^d of the spacecrafts - i.e. both terms of equation (16) have 

negative values. The same restriction applies to equation (16) as to equation 

(10), namely that a^ is not a modification of the gravitational acceleration 

due to the Sun. The average radial velocities for Pioneer 10 at distances 
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Table II: Predicted values for the Pioneer Anomaly acceleration. 



Distance (AU) 


aA X 10^" ms- 


-2 


Pioneer 11 




22 


11.7 


27 


10.1 


32 


9.2 


Pioneer 10 




40 


8.5 


50 


7.9 


60 


7.6 


70 


7.4 


Average 


9.0 



from 40 AU to 70 AU and for Pioneer 11 at distances from 22 AU to 32 AU 
were, respectively, 12.6 kms"""^ and 11.8 kms^^ ^ . Substituting these values 
in equation (16) and using a value for H of 71 kms^^ Mpc~^ and Mq — 
2 X 10'^° Kg generates values for aA as shown in Table II. 

All the predicted figures for aA lie within the tolerance limits reported 
for the Pioneer Anomaly, 8.74± 1.33 x lO^^'^ ms^^, except for those for Pio- 
neer 11 at distances less than 27 AU. However, the most significant feature 
of equation (16) is, perhaps, that it predicts a variable value for a a rather 
than the reported constant value. 



5. Discussion 

An analytical framework based on simple assumptions about four dimen- 
sional space-time being closed, isotopic and expanding, and embedded in a 
fifth large-scale dimension, r, representing the radius of curvature of space- 
time, was originally developed by the authors as a tool to study the orbital 
dynamics of galaxies [3] . The analysis presented in this paper implies that 
this framework may have wider applications. In both cases, the principles 
of Newtonian dynamics and gravity are applied without amendment. 

For the predicted additional acceleration [a a) to give rise to an observ- 
able effect on a body (i.e. observable within four dimensional space-time), 
it must be moving within a system which permits transfer of kinetic energy 
to or away from it, such as through planetary flybys. Moreover, in the case 
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of flybys the net effect of a^ over the combination of inbound and outbound 
trajectories will only be measurable to the extent that these trajectories are 
asymmetrical with respect to their energy transfer characteristics. 

Several parameters, such as the ecliptic angle of incidence of the trajec- 
tory and its radial velocity towards (or away from) the planet, determine 
this asymmetry; however, as Anderson et al.^ have shown, the equatorial 
latitude appears to be a good surrogate for this measure of asymmetry, at 
least in the case of the six Earth flyby trajectories considered. The central 
role of radial velocity in determining ua for Earth flybys (equation 10) sug- 
gests why no such velocity discrepancies have been detected for satellites 
in bound orbits. 

The aggregate force, F, due to the combination of gravitational and 
the additional acceleration term acting on a mass, m, having radial motion 
relative to the gravitating body is, using equation 11, given by equation 
17, the form of which has an intriguing similarity to the Lorcntz force of 
electromagnetism. 

F = m (gs - - A A (17) 

The equations governing aA have been derived from first principles and 
without the introduction of free parameters or new constants. That they 
appear potentially to provide a single explanation for both Earth flyby and 
Pioneer anomalies is an added attraction and justifies their further testing 
through application in detailed dynamic models of the solar system. 

6. Addendum 

6.1. Introduction 

Although initially promising agreement was achieved between the analyt- 
ical framework proposed in this paper and the reported fiyby anomalies, 
subsequent detailed modelling undertaken with the assistance of Anderson, 
suggested the agreement to be much poorer than expected, especially when 
specific aspects of the Doppler tracking data were considered. Furthermore, 
the issue of an unidentified energy transfer mechanism (required, if the 
anomalous acceleration is to be real) remains open, whether energy trans- 
fer is correlated to the geocentric latitudes of the spacecrafts inbound and 
outbound velocity vectors (Anderson et al), so implying that the Earth's 
spin energy is the source/sink for energy transferred, or whether the Earth's 
orbital energy around the Sun is proposed as the source/sink, as in this pa- 
per. Evaluation of these results and the issue of energy transfer has led to 



April 12, 2010 0:34 WSPC/Trim Size: 9in x 6in for Proceedings gerrard3v2 



11 

the conclusion that the acceleration, a a, derived from consideration of five 
dimensional dynamics could, alternatively, be interpreted as virtual rather 
than real. That is, interpreted as an acceleration that affects the orienta- 
tion of spacecrafts momentum vector but not its energy, so influencing its 
direction of motion relative to the Earth, rather than its speed. It is the 
purpose of this addendum to show how this alternative interpretation of aA 
may overcome the principal short-comings so far found in the ability of five 
dimensional dynamics to model accurately the reported flyby anomalies. 



6.2. Interpreting the additional acceleration, a a, as virtual 

In the case of Earth flybys, the additional acceleration, a a, experienced by 
the spacecraft is directed towards the Earth, in accordance with equation 
(10), even if it is virtual in nature. However, being virtual it can be treated 
as analogous to a centripetal acceleration affecting the hyperbolic orbit of 
the spacecraft about the Earth. Conservation of energy then requires that 
the effect of aA on the spacecraft will be such as to induce a change in direc- 
tion of its motion (effectively a perturbation of the spacecraft's momentum 
vector) such that the resultant induced small increase in centrifugal compo- 
nent of acceleration of the spacecraft about the Earth (defined as Ucf) will 
be equal and opposite to a a- This principle is illustrated mathematically 
(and in its simplest form) in equation (18) and diagrammatically in Figure 
4. Clearly, however, to the extent that the trajectory of the spacecraft, 
shown as mass m, departs under the influence of ua from that expected in 
four dimensional dynamics, its time-dependent profile of kinetic energy and 
gravitational potential energy inter-change will itself be perturbed and, to 
that extent, a a would be expected indirectly to affect the time-dependent 
kinetic energy profile of the spacecraft as it approaches/recedes from Earth, 
as a real although second order effect. 



aA + dttcf = (18) 

Figure 4 considers a point on the trajectory of a spacecraft (mass m) in 
motion towards a flyby of the Earth (mass M) , and shows the separation 
of the spacecraft's velocity, v, into two orthogonal components: one, Vradi 
directed radially towards the Earth and the other, Vtam normal to the 
radial vector, s. The centrifugal component of acceleration acting on the 
spacecraft with respect to the Earth is given by equation (19), where s is 
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Figure 4. Tangential velocity perturbation, dvtan, induced by a^- 



the distance between spacecraft and Earth: 

acf = "^ (19) 

s 

The proposed effect of qa is to induce an increase in tangential compo- 
nent of velocity, dvtam of the spacecraft. The gain in tangential component 
of velocity, dvtam induced by a a will itself be offset by a reduction in radial 
component of velocity of the spacecraft with respect to the Earth, dvrad, for 
energy to be conserved under the action of a a, in accordance with equation 
(20), for vtan » dvtan and Vrad » dvrad- The relationship between an in- 
cremental increase in induced centrifugal component of acceleration, dacf , 
and the associated induced incremental increase in tangential component 
of velocity, dvtam is given by equation (21): 

VtandVtan + VraddVrad = (20) 

[acf + dacf) = (21) 

From equations (10), (18), (19), (20), (21) and substituting the grav- 
itational acceleration vector, (/s, by its component form of g^ = G.M/s'^, 
for Vtan » dvtan and for Vrad » dvrad > 0, it cau be shown that the 
corresponding incremental decrease in the spacecraft radial component of 
velocity, dvrad, is given by the equation: 

dVrad^-^ (22) 

2cs 
Care is needed in applying equations (22) because it is valid only for 
constant spacecraft kinetic energy, as it ignores the time-dependent effect of 
the Earths gravitation acceleration on v, Vrad and Vtam in accordance with 



April 12, 2010 0:34 WSPC/Trim Size: 9in x 6in for Proceedings gerrard3v2 



13 

four dimensional dynamics. Its use is, therefore, strictly limited to being an 
equation of instantaneous perturbation between four dimensional and five 
dimensional dynamics. This limitation on equation (22) is especially rele- 
vant in the region at or about the flyby itself when, from equation (10), it 
is apparent that as Vrad tends to zero, so does ua which means that pertur- 
bations to the tangential and radial velocity components of the spacecraft, 
d'>^tan and dvrad wiU also tend to zero at the spacecraft perigee. Also, in 
applying equation (22) it should be noted that we are concerned only with 
positive values for dvtan and so only with negative values of dv^ad- More- 
over, the relative contributions of Vrad and vtan towards the spacecrafts 
velocity, v, vary continuously along the length the spacecrafts inbound and 
outbound trajectories according to the distance s and angle 6 subtended 
by the velocity vector v with respect the radial vector. For these reasons 
numerical integration is required to model accurately the overall effect of 
GA on spacecraft trajectory, which is beyond the scope of this paper. How- 
ever, from equation (22) it is possible to see the instantaneous perturbing 
effect that ua as a virtual acceleration would have on the spacecrafts tra- 
jectory. The spacecrafts inbound trajectory will be slightly skewed away 
from Earth relative to a trajectory derived by reference solely to four di- 
mensional dynamics (due to a slightly accelerated transition to the purely 
tangential velocity condition at the perigee), leading to a slightly higher 
altitude perigee than expected. Conversely, the outbound trajectory will 
be slightly skewed towards the Earth relative to a trajectory derived by ref- 
erence solely to four dimensional dynamics. This implies that the perturb- 
ing effect of ga on spacecraft trajectory would, therefore, be asymmetrical 
about the perigee. 

Lastly, it needs to be stressed that in this analysis the perturbation of 
tangential velocity, dvtam induced by a a has been assumed to lie in the 
direction of the tangential velocity vector itself, Vtan ■ This need not neces- 
sarily be so, as equation (18) could, in principle, equally well be satisfied by 
a component of perturbed tangential velocity lying anywhere in the plane 
normal to the radial vector, s, between spacecraft and Earth. Hence, in 
the limiting case, equation (18) could even, in principle, be satisfied by an 
induced gradual gyration of the spacecraft about the radial vector. If this 
were the case, it may not be directly apparent from the Doppler tracking 
data which measure spacecraft range and radial velocity only. However, 
it seems reasonable to proceed on the basis of there being a preferred di- 
rection for dvtam not least because the radial motion of spacecraft with 
respect to the Sun's gravitational field will itself give rise to a small com- 
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ponent of a Sun-directed additional and virtual acceleration, a a, which 
defines a preferred direction for dvfan relative to the Earth. Finally, it is 
important to note that equation (18) is not the only possible interpretation 
of how a virtual acceleration, a^, could induce changes in the direction of 
the spacecraft velocity vector. It is chosen here for illustration because it 
is mathematically the simplest. 



6.3. Discussion 

The above analysis illustrates how a a, although virtual in nature, can give 
rise to small although systematic and asymmetric perturbations of the in- 
bound and outbound trajectories of a spacecraft during Earth flybys. Since 
the inference of anomalous changes in spacecraft energy during flybys re- 
quires the observed radio Doppler tracking data to be fitted to equations of 
four dimensional dynamics, it is apparent that such unanticipated perturba- 
tions may potentially give rise to apparent anomalies. There are potentially 
two important sets of unknowns in the spacecraft trajectories, as defined by 
the Doppler tracking data: the first is the tangential component of velocity, 
Vtan which is not directly determined; and the second derives from any ab- 
sences of reliable Doppler data at or about the perigee. The potential loss 
of Doppler data during the perigee may derive from a combination of the 
latitude of spacecraft relative to available tracking stations and difficulties 
in data quality that can arise from tropospheric refraction. For example, 
the gaps in Doppler data for both NEAR(1998)i and Rosetta(2005)'^ fly- 
bys are over 1 hour in duration (not necessarily symmetrically centred on 
the perigee). Both sets of unknowns are deduced from fits of the avail- 
able Doppler data to the equations of four dimensional dynamics. Since 
the essence of the proposed additional virtual acceleration, a a, is to induce 
both a perturbation to the tangential component of spacecraft velocity, 
dvtam and a slightly higher altitude perigee than expected, there is clearly 
potential for it to give rise to apparent anomalies in data fitting. The task 
of using Doppler tracking data to infer the extent of trajectory perturbation 
introduced by qa niay be made harder than expected by the functional form 
of equation (22). This equation shows that dvrad is proportional to GM/s 
and so implies that the distance-related dv^ad profile of the spacecraft due 
to five dimensional effects, should be similar to the distance-related profile 
of spacecraft gravitational potential energy under four dimensional dynam- 
ics, which is also proportional to GM/s. The point can be illustrated by 
considering sections of pre and post encounter trajectories where the ap- 
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proximation v « Vrad is valid. In these cases, the incremental change in 
spacecraft velocity (Av) at a point along the trajectory can be approxi- 
mated by the addition of two components: one due to four dimensional 
dynamics (Avid)] and the other due to five dimensional effects (Av^d)'^ 
the latter being taken from equation (22). This approximation is shown in 
equation (23) and is rearranged in equation (24)'', for s >> As, to highlight 
the way in which five dimensional effects would show-up as a perturbation 
(cp) to the four dimensional case. 

Aw « Av4d + Av5d ~ Av4,d (1 + fp) (23) 

GMAs / w \ 

In fitting the Doppler data to four dimensional equations of motion, a time 
dependent perturbation introduced by the term v/2c may be absorbed by 
the imputed changing values of Av (and so v) as an underestimate of the 
implied value of the tangential component of velocity of the spacecraft, 
which is not independently measurable. Depending on the step length, 
As, used in the numerical integration, the perturbation could fall within 
the error bars of the methodology where the approximation v « Vrad is 
valid. For example, applying equation (24) for a spacecraft at a distance 
from Earth, s, of 10^ km, an integration step length. As, of ^ 10'^ km or 
less (equivalent to 100 seconds or less, at a velocity of 10 kms~^) would be 
consistent with a fit to four dimensional dynamics within an equivalent noise 
level of less than ~0.1mins^^ in spacecraft velocity. Consequently, both 
inbound and outbound limbs of the spacecrafts trajectory could provide 
good fits to four dimensional dynamics and yet be inconsistent with each 
other. This inconsistency could show up, inter alia, as an un-modellcd 
change in spacecrafts energy. Such an effect would be consistent with the 
anomaly becoming immediately apparent, within the fitted trajectories, 
from the point at which the spacecraft emerges from its encounter with the 
Earth and radio tracking contact is resumed'. 

Moreover, the rate of decline of qa with spacecraft distance from Earth 
evident from equation (10), and the inverse dependency of Avrad on space- 
craft distance from Earth evident in equation (22) would imply a strong 



Av5£j Si dv^ad (* ~ ^*) ~ d,Vj.ad (s) from equation (22) 
''In equation (24), the sign will be reversed for outbound trajectories 

A similar effect will be expected in fitting trajectories to 4-dimensional dynamics if aA 
were interpreted as real and Av^d were taken from equation (13) rather than (22) 
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distance dependency for difficulties in achieving consistent fits to four di- 
mensional dynamics, in cases where the proposed five dimensional pertur- 
bation effects arise. The consequent reducing potential for reported anoma- 
lies in cases of higher altitude flybys is consistent with the conclusions of 
Anderson and Nieto [7]. 

In particular, equations (23) and (24) can be rearranged to estimate the 
altitude, Humit, below which the perturbation Aw5d will exceed a specified 
sensitivity level (e.g. O.Imms^^) as follows: 

^H™*«(^^^^J -REartH (25) 

For a spacecraft velocity ^ lOkms^^ and integration step length of Is, 
Hiimit is then ^ 2000 km. For a given satellite Humit will depend on its 
detailed trajectory. However equation (25) does indicate the altitude above 
which anomalies would not be detectable, and is in general agreement with 
observation to date [7]. 

It can be further noted that whereas the interpretation of aA as real, 
and acting in a radial direction towards Earth, may present difficulties in 
interpreting anomalous changes in spacecraft angular momentum during 
the fiyby, the interpretation of a a as virtual may not, as it concerns the 
induction of perturbations to the spacecraft velocity vector but not its 
energy. 

The interpretation of aA as virtual is not specific to Earth flybys, so 
some comment can also be made about the Pioneer Anomaly and equation 
(16). This equation has two components, the first {gsVrad/c) derived from 
five dimensional dynamics, and the second {—cH) derived from considera- 
tion of the background expansion of a five dimensional universe. Whereas 
the first component has a clear virtual analogue inducing changes in space- 
craft momentum vector but not in its energy (so causing the instantaneous 
tangential component of Pioneer spacecraft velocity to increase at the ex- 
pense of the radial component), the second component of equation (16) 
does not. Considering only the first term, therefore, and for a spacecraft 
on a hyperbolic orbit travelling away from the Sun, it is possible to use 
equation (22) to derive a time derivative of the perturbation to the radial 
velocity, dv^ad, of the spacecraft that is induced by the virtual acceleration, 
a A, relative to that expected under four-dimensional dynamics. The accel- 
eration, a A, is in this case directed towards the Sun. The time derivative 
is given by equation (26) which is simply half the effect of that implied if 
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a^ were treated as real (i.e. half the first term of equation (16)). 

dt ~ 2c~ 

The values of the Sun-directed acceleration computed from equation (26) 
for Pioneers 10 and 11, lie in the range 0.3 to 2.4 x 10^^" ms^^ which is well 
below the reported values implying that this effect alone cannot account 
for the anomaly [8,9]. 

6.4. Conclusions 

This Addendum has been included to address the possibility of the ad- 
ditional acceleration, a a, derived from five dimensional dynamics, being 
interpreted as virtual rather than real. The analytical framework used is 
otherwise unchanged from that in the body of this paper. The next step, as 
before, is to carry out numerical integration of spacecraft equations of mo- 
tion modified for the perturbing effects of aA and to reconcile these results 
with trajectory fits to available Doppler tracking data, to see whether these 
perturbations are consistent with the reported flyby and Pioneer anomalies. 
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